Non-linear tx-model for long range disorder and quantum chaos. 



(N 
O 

o 

(N 

> 
O 



S3 
S3 



V. R. Kogan 1 ' 2 and K. B. Efetov 1 * 2 
1 Theoretische Physik III, 
Ruhr-Universitat Bochum, 44780 Bochum, Germany 
2 L. D. Landau Institute for Theoretical Physics, 117940 Moscow, Russia 

(February 1, 2008) 

We suggest a new scheme of derivation of a non-linear ballistic cr-model for a long range disor- 
der and quantum billiards. The derivation is based on writing equations for quasiclassical Green 
functions for a fixed long range potential and exact represention of their solutions in terms of func- 
tional integrals over supermatrices Q with the constraint Q 2 — 1. Averaging over the long range 
disorder or energ y we are able to write a ballistic cr-model for all distances exceeding the electron 
wavelength (Eq. 3.13). Neither singling out slow modes nor a saddle-point approximation are used 



in the derivation. Carrying out a course graining procedure that allows us to get rid off scales in 
the Lapunov region we come to a reduced cr-model containing a conventional collision term. For 
quantum billiards, we demonstrate that, at not very low frequencies, one can reduce the cr-model 
to a one-dimensional a-model on periodic orbits. Solving the latter model, first approximately and 
then exactly, we resolve the problem of repetitions. 

PACS: 05.45.Mt, 73.23.-b, 73.23.Ad 
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I. INTRODUCTION 

The cr-model approach to disordered systems, first 
written within the replica tricktlu, proved to be a power- 
ful method of calculations when formulated in the su- 
persymmetric formEI. This method allows to describe 
the electron motion at large distances assuming that at 
shorter distances the motion is diffusive. For study of 
such phenomena as localization, level statistics in a sys- 
tem of a size L much exceeding the elastic mean free 
path, etc., the information obtained from the cr-model is 
sufficient. Although the cr-model is not valid at distances 
smaller than the mean free path I, the information about 
the motion at distances below I is not very interesting for 
these phenomena. 

Success in nano-fabrication made possible producing 
and studying clean systems with the size smaller than the 
elastic mean free path I (for a review, see. e.g.cl). In such 
small systems called now quantum dots electrons move 
ballistically being scattered mainly the walls. There are 
many interesting questions about transport in the quan- 
tum dots and related systems that cannot be answered 
using the picture of the diffusive motion. In order to de- 
scribe the ballistic motion one has to go beyond the di- 
agrammatic and field theoretical methods developed for 
disordered systems. 

Another motivation to study the ballistic motion origi- 
nates from the field called now quantum chaos. The sub- 
ject of research in the quantum chaos is to understand 
the quantum behavior within models that are chaotic 
in the classical limit. There are iflpftny books and re- 
views related to this field (see, e.g.BEm). The most pop- 
ular analytical tool for studying the quantum chaos is 
the Gutzwiller trace formulatl that reduces calculation 
of the density of states to a sum over periodic orbits. 



This method (complemented by different approximation 
schemes) allows one to study very well the limit of not 
very long times t when the motion along a periodic orbit 
is well defined. At the same time, calculations with the 
trace formulae become very difficult in the limit t — > oo 
when one expects an universal behavior described by the 
Wigner-Dyson statistics^. 

Both the experimental and theoretical interest to in- 
vestigations of the ballistic motion resulted in several at- 
tempts to construct a generalization of the supersymmet- 
ric cr-model to distances smaller than the mean free path 
I due to scattering on impurities or the wallsJn the sys- 
tem. Muzykantskii and Khmclnitskii (MK)li3decoupled 
as usual the i/> 4 interaction in the effective Lagrangian 
by gaussian integration over a supermatrix Q but did 
not use after that a saddle-point approximation. Instead, 
they derived a quasi-classical equation for effectiveXlreen 
functions g analogous to the Eilenberger equationU well 
known in the superconductivity theory. Using an analogy 
of this equation with an equation for motion of a mag- 
netic moment in an external magnetic field MK noticed 
that the equation was a minimum of a functional <I> con- 
taining a Wess-Zumino-Novikov-Witten (WZNW) term. 
So, they replaced the solution of the semiclassical equa- 
tions by a functional integral containing the functional $, 
which allowed to average over the supermatrix Q. This 
could be done provided the equation of the Green func- 
tions corresponded to a deep minimum of the functional 
$, such that fluctuations near the minimum could be ne- 
glected. Although the authors of Refill conjectured that 
their field theory could be applicable even in the limit of 
a vanishing disorder, they did not confirm this point of 
view by any calculations. 

A more traditional way of derivation was used by An- 
dreev et af-Hl who tried to derive the a-model for a bal- 
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listic quantum billiard. Instead of averaging over dis- 
order they averaged over the energy. After decoupling 
the ^ 4 term by integration over the supermatrix Q they 
used the saddle-point approximation, which is equivalent 
to the self-consistent Born approximation (SCBA). This 
had to fix the eigenvalues of Q such that one could put 
Q 2 = 1 and what remained to do was to expand the ac- 
tion in gradients of Q and in the frequency u). Shortly 
after it became clear that the saddle-point did not fix the 
eigenvalues of the supermatrix Q and modes that were 
usually massive for disordered system became massless 
in the ballis tic , li mit. This problem was discussed in the 
publicationst^rEJ. It is relevant to notice that the ex- 
istence of the additional massless modes is not just a 
consequence of the bad saddle-point approximation. As 
discussed in Ref .113 the same problem is encountered when 
deriving the cr-model with the help of the so-called "color- 
flavor" transformation. Within this approach, although 
one does not need to use any saddle-point approximation, 
the expansion in gradients still remains to be performed. 
However, there is no parameter that would allow one to 
take into .account the lowest gradients only. Zirnbauer 
suggestcdtij to perform an additional averaging over en- 
sembles in order to suppress short range fluctuations. 

The saddle-point approximation can actually be useful 
if one considers a system with a long range disordered. 
In this case, the single-particle mean free path I can be 
much smaller than the transport mean free path l tr . At 
distances exceeding the length I the saddle-point approx- 
imation and the expansion in gradients of Q can be used 
and one comes to a ballistic cr-model that reduces to the 
diffusion one only at distances exceeding l tr . Thus, in 
the interval between I and lt r one can obtain the ballistic 
cr-model in a reliable way (see also a subsequent discus- 
sion in Ref.t3). However, this does not solve the problem 
completely because a reasonable semi-classics should be 
applicable at all distances exceeding the wavelength \p. 

A ballistic cr-model should describe low lying excita- 
tions that exist for any long range disorder. At the same 
time, the conventional saddle-point approximation, being 
equivalent to the SCBA, can be good for a short range 
disorder only and, hence, may not be used for deriva- 
tion of a cr-model for a long range disorder and quan- 
tum chaos at arbitrary distances. The same is true for 
the decoupling of the ip 4 term. The integration over the 
supermatrix Q is usually used after singling out slowly 
varying pairs ipip. However, if the random potential is 
very long ranged or one averages over the energy, one 
has slowly varying pairs from the beginning and there 
is no necessity of integration over the supermatrix Q 
instead of integration over the initial random potential 
U (r). The sarne is true when applying the color-flavor 
transformationEj. The replacement of an integration over 
u (r) by an integration over a supermatrix Z (r) does not 
seem to correspond to physical processes and is a purely 
(although exact) mathematical transformation. Thus, 
the correct scheme of the derivation of a field theory de- 



scribing the low lying excitations should not be based on 
the Hubbard-Stratonovich or color-flavor decoupling and 
the saddle-point approximation. 

In this paper, we present a derivation of a ballistic 
cr-model using neither the Hubbard-Stratonovich decou- 
pling with a supermatrix Q nor the saddle-point approxi- 
mation determining the eigenvalues of Q. As in Ref.E2l we 
derive quasiclassical equations for Green functions but we 
write them without making the Hubbard-Stratonovich 
transformation. This is justified because we use a long 
range potential. Only if a short range potential is added 
we must single out slow pairs in corresponding term in 
the Lagrangian and decouple it by the integration over 
supermatrices. The crucial step of the derivation is an 
exact representation of the solution of the quasiclassical 
equation in terms of a functional integral over 8x8 su- 
permatrices Q n (r) with the constraint Q\ (r) = 1, where 
r is the coordinate and n = p F / \pp\ is the normalized 
vector on the Fermi-surface. An effective action $ u [<?„] 
entering, the functional integral is similar the one written 
in Ref.li3. We show that supersymmetric properties of the 
matrix Q n (r) make the representation exact, which was 
not noticed in the MK variational approach. Moreover, 
the solution written for an arbitrary long range potential 
u (r) is applicable even for non-averaged quantities. 

Averaging over the random potential leads to an ef- 
fective action $ [g n ] that has a form different from those 
discussed previously. Analyzing properties of the new 
ballistic non-linear cr-model with the action $ [g n ] we 
demonstrate that a new length li, = vpr^ introduced 
by Aleiner and LarkinliU, where vp is the Fermi velocity 
and tl is the inverse Lapunov exponent, determines dif- 
ferent regimes. The importance of this length was also 
discussed recently in Ref £3. Integrating over variations of 
the supermatrix Q n (r) at distances smaller than lp we 
come to another form of the ballistic cr-model containing 
the conventional collision term. 

We show that without an internal disorder the calcula- 
tion of the functional integral can be reduced to study of 
the cr-model for periodic orbits. Only the presence of a 
regularizer analogous to the one introduced in Ref.EZl.may 
mix the periodic orbits. The problem of repetitionalij is 
discussed and we are able to demonstrate that the con- 
tradiction between the references^ ano£3 is rather a cop. 
sequence of an unjustified approximation used in Ref.Eij 
than a deficiency of the cr-model. 

The paper is organized as follows: In Chapter II, we 
express correlation functions of interest in terms of func- 
tional integrals over supervectors and write equations for 
generalized Green functions. In Chapter III, we represent 
the solution of the quasiclassical equations in terms of 
functional integrals over supermatrices and average over 
disorder, thus obtaining a ballistic cr-model applicable at 
all distances exceeding the wavelength. In Chapter IV, 
we integrate over a Lapunov region and derive a reduced 
ballistic cr-model containing a collision term. In Chapter 
V, we show how one can derive equations for correla- 
tion functions. In Chapter VI, we show how calculations 
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within the ballistic cr-model can be reduced to calcula- 
tions for periodic orbits. We explain how the so called 
"repetition problem" can be resolved. Chapter VII is 
devoted to a discussion of the results obtained. The Ap- 
pendix contains a derivation of the boundary conditions. 



II. FORMULATION OF THE PROBLEM. 
QUASICLASSICAL APPROXIMATION. 

The aim of the present paper is to find a convenient 
representation that would allow us to consider an elec- 
tron motion in a smooth potential at large times or low 
frequencies. Of course, with the formalism presented one 
can consider wave scattering in microwave cavities and 
other interesting problems but, to simplify notations, we 
will use the condensed matter language. n 

We want to extend the supersymmetry methods devel- 
oped for disordered systems to distances smaller than the 
mean free path. Actually, the only assumption we will 
use in the derivation is that all physical quantities vary at 
distances exceeding the Fermi wavelength Xp = 2-Kpp 1 ' . 
Although this assumption is much less restrictive than 
those used in Ref.cl it allows to simplify essentially the 
consideration. Averaging over the energy, which is the 
standard procedure for quantum chaos, can be consid- 
ered as the limiting case for an infinite range random 
potential. 

The method developed in this paper is applicable for 
calculation of gauge invariant quantities like density- 
density or level-level correlation functions. Such quan- 
tities as average one-particle Green functions at different 
points will not be considered here. We choose the Hamil- 
tonian H of the system in the standard form 



H = Hq + u (r) + u s (r) , 



(2.1) 



H 



-V 2 /2m-e F 



where u (r) is a long range potential, which is of the main 
interest now, and u s (r) is a short range impurity poten- 
tial. The latter is added in order to make the models 
somewhat more general. The presence of the short range 
potential will help to understand better the procedure we 
will use. However, nothing is assumed about the strength 
of u s (r) and it can be safely put to zero in all formulae 
written belgw. 

As usualU, one can express correlation functions of in- 
terest in terms of a functional integral over 8-component 
supervectors V (r) with an effective Lagrangian L 



L bP] 



iip(r) (H + u(r))ip(r) 



i (lu + id) 



V (r) Atp (r) 
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(V(r)V(r))V, (2.2) 



Ho — H - e + — 

where e is the energy at which the physical quantities are 
considered and w is the fre que ncy. 

The Lagrangian L, Eq. ( |2.2| ), is written after the aver- 
aging over the short range potential u s (r) implying the 
standard gaussian correlations of the type 



(u s (r) u s (r')) 



1 



27WT. 



■J(r-r') 



(2.3) 



It is relevant to emphasize that averaging over the long 
range potential u (r) has not been performed. 

The correlation functions we are interested in can be 
obtained adding proper source terms in the Lagrangian 
L [tp] . An important class of the correlation functions can 
be obtained writing the Lagrangian L a [tp] including the 
sources in the form 



L a [ip]=L[ip]+i / V (r) a (r) j> (r) dr 



(2.4) 



where a (r) is a matrix depending on coordinates. Its 
explicit form depends on what type of the correlation 
function is calculated. 

The level-level correlation function R (u>) 
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2tt 2 uj V 2 V 2 



(Re / (n (e — uS) — n (e)) 



xGt,, (r,r) (Gf (r',r') - (r',r')) dvdv'de), (2.5) 

where n (e) is the Fermi distribution, contains the prod- 
uct Gf^Gf. For calculation of this product one should 
choose a in the form 



&! 
—&2 



a la = ^f{l-k) (2.6) 



where both k and T3 denoting different blocks have the 
form 

' 1 
-1 / 

Below we use the same notations as in the book!. 

The product (Gf _ M Gf ) is trivial, whereas the first 
product in Eq. (2.5) can be written as 



(G?_ u (t,t)G? (r',r')) 



4 / VI (r) VI (r) VI (O VI (O exp (-L a [V]) Dxj) 



and we reduce the function R (u>) to the form 
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(2.7) 
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lira Rc 



2(tTVV) «1=«2=0 



d 2 



da\da2 



exp(-L a [ip})Dtl> 



where 



For computation of the density-density correlation func- 
tion one has to calculate the averages of the type 



J (r) = ia (r) + 



M(r) 
2t„ 



(2.12) 



y 00 (r 1 ,r 2 ;^) = 2(G^(r 2 ,r 1 )Gf (r 1; r 2 )) (2.8) and L, 7 [^ is obtained from L a [V], Eqs. @, @, after 



This product can be obtained from the following source 
term 



where 



a (r) 



«o (r) 
-a (r) 



&i<5(r — ri) 

a 2 5(r — r 2 ) 



(2.9) 



with the same expression for 6t\,2 as in Eq. (2.6). Then, 
we have for F 00 (ri, r 2 , cj) 



F 00 (r 1 ,r 2 ; W ) = 2 f e xp 



«1— «2— 



(2.10) 



In principle we can proceed calculating in Eqs. (E 



2.10) in the standard waya by singling out slowly vary- 
ing pairs tpij) in the term ip 4 in the Lagrangian L, Eqs. 
(2.2, 2.4), and decoupling the products of these terms 



by Gaussian integration over 8x8 supermatrices M. 
A word of caution should be said at this point. The 
separation into the products of slowly varying pairs ipip 
makes a sense only for distances exceeding the range of 
the random potential. This means that one may not this 
approximation for distances smaller than the potential 
range. Of course, the same is true when averaging over 
the spectrum. In the latter case the range of the random 
potential is just the system size. Thereforej-the previous 
derivations where this separation was usedt3 can hardly 
be justified. The same problem apparently arises when 
doing the color-flavor transformation^. Although the 
transformation is formally exact, there is no reason for 
neglecting higher gradients when making the expansions 
in gradients. 

In order to avoid the problem we keep the long range 
potential u (r) in Eq. (pO,2.4) as it stands and do not av- 
erage over it. This will be done later. At the same time, 
the decoupling of the term corresponding to the short 
range potential by integration over the supermatrix AI 
can be safely done (We repeat again that the presence of 
the short range potential is not crucial for our derivation 
and its strength can be put zero). 

After the decoupling for the short range potential the 
calculation of the correlation functions is reduced to the 
computation of an effective partition function Z\ [J] 



Zi [J] 



J exp(- 



Di\) 



(2.11) 



the r eplac ement of the parameter a by J according to 
Eq. ( 2.12 ) and neglecting the quartic term. The func- 
tion J satisfies the standard symmetry relation J = J, 
where J = CJ T C T .r-,The bar means the usual "charge 
conjugation" of Refs.0, the matrix C is defined as follows 



C = A< 



ci 
c 2 



ci = 



-1 

1 



C2 = 



1 

1 



The correlation functions for a given long range poten- 
tial u (r) can be calculated by differentiating in a.x,2 the 
following integral 

Z u = J Zi [J] exp \-j£r J StrM 2 (r) dr^j DM (2.13) 

Averaged correlation functions can be obtained from the 
quantity Z = (Z) u , where (. . .) u means averaging over 
u (r) . Of cou rse, we could immediately average over u (r) 
in Eq. ( 2.13| ) but we want to avoid the standard scheme. 
The approximations that worked so well for short range 
potential are not applicable to the long range one. In 
particular, not only the separation of slow modes is not 
justified but also the saddle-point approximation for the 
integral over supermatrices is no longer good. It is clear 
that the existence of low lying excitations like diffusons 
and cooperons is more general than the SCBA and one 
should try to avoid the latter. 

In this paper we follow the method of quasiclas- 
sical Green functions first introduced for study of 
superconductivity^. This method was used recently by 
Muzykantskii and Khmelnitskiilij for study of the ballis- 
tic transport -i-Our calculations are partially equivalent to 
those by MKE3 but there are essential differences. First, 
we keep the long range potential u (r) fixed and aver- 
age over it at the later stage of the derivation. Second, 
which is the most crucial step, we show how to write the 
solution of the quasiclassical equations exactly. The pos- 
sibility of writing an exact solution of the equations for 
quasiclassical Green functions in terms of functional in- 
tegrals is a consequence of the supersymmetry and has 
not been realized before. 

As in MKLL3, we consider the logarithmic derivative of 
the partition function Z\ [J (r)] 



(51nZi[J(r)] 
5J al3 (r) 



(2.14) 



where < . . . >^ is th e average with the functional Lj [ip] , 
Eqs. (ET|, 121 HI) and 
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Introducing the matrix function G (r, r') as 

G(r,r') = 2(V-(r)Vi(r'))^ (2.15) 
we can write for this function the following equation 



~ / „ , uj + iS 

H Qr + u (r) + A—— + U (r) 



G(r,r') = iJ(r-r') 
(2.16) 



where the subscript r of H$ r means that the operator 
acts on r. 

Conjugating Eq.( 2.16j ) and using the property 



G(r, r') = GG T (r', r)G T = G (r, r') 



(2.17) 



we obtain another equation for the matrix G(r, r') with 
the operator Ho r > acting on its second variable 



G(r;r') 



= z<S(r-r') 
(2.18) 



Until now no approximations have been done and Eqs. 
(|2.16, 2.18|) are exact. Now we can use the assumption 



that the potential u (r) changes slowly on the wavelength 
Xf . If the mean free path for the scattering on the ran- 
dom potential exceeds the Green function varies as 
a function of r — r' at distances of the order of Xf but, 
at the same time, is a slow function of (r + r') /2. The 
Fourier transform G p ((r + r') /2) of G (r, r') respective 
to r — r' has a sharp maximum near the Fermi surfa ce. In 
order to c ancel large terms we subtract Eq. ( [2.18 ) from 
Eq. ( |2.16 ). Using the assumption that the potential u (r) 
is smooth and expanding it in gradients we obtain in the 
lowest order 



--^-^+zV R u(R) — 
m op 



G P (R) 



(2.19) 



id 



[A,Gp(R)]+i[J(R),G p (R)] = 0, 



where R = (r + r') /2 and 



stands for the commutator. 
When deriving Eq. (2.19), not only the potential u (r) 



but also the function J (r) was assumed to be smooth. 
The supermatrix M (r) is smooth by the construction. 
As concerns a, this is not always so, as is seen from Eq. 
(2.8-2.10). However, we can slightly smear the coordi- 
nates in the definition of the correlation functions and 
obtain after this procedure a smooth function a. 

The dependence of the Green function G p (R) on |p| is 
more sharp than on other variables. I n orde r to avoid this 
sharp dependence we integrate Eq. ( 2.19| ) over |p|. Of 



course, this procedure makes a sense for very large sam- 
ples when the level discreteness can be neglected. How- 
ever, this procedure can also be performed in finite sam- 
ples provided it is complemented by an averaging over 
the energy. 

The most interesting contribution in the integral over 
|p| comes from the vicinity of the Fermi-surface. A con- 
tribution given by momenta considerably different from 
Pf is prop ortional to the unity matrix an drops out from 
Eq. (|2~19|) . 

Introducing the function g n (r) 



7n (r) = i [ G pn (r)d£, £ I- ^ 



(2.20) 



where n is a unite vector pointing a direction on the 
Fermi surface, we obtain the final quasiclassical equation 

{v F nV -p^V r u(r)d n ) gil {r) (2.21) 

i(u + iS) 



where 



-[A,£ n (r)]-[J,<7 n ]=0 



9n = V r 



d 

V n = -[n x [nx— ]] 
on 

The function g n (r) is self-conjugate 
g n (r) = Cg T n (r) C T — g n (r) , 



(2.22) 



which follows from Eqs. fl2.17| ), fl2.20| ). 

Eq. ( 2.21 ) should be complemented by a boundary 
condition at the surface of the sample. This boundary 
condition is derived in Appendix A. Considering a closed 
sample we assume that the current across the border is 
equal to zero. This leads the boundary condition at the 
surface 



?nj_ ( r )\ SU rface — 9-n± ( r ) 



surface 



(2.23) 



where nj_ is the component of the vector n perpendicular 
to the surface. 

The correlation functions considered here can easily be 
expressed in terms of the quasiclassical Green functions 
(? n (r) . The f unctio nal derivative of the partition func- 
tions Zi, Eq. ( |2.14 ), can be written through the function 
g n (r) as follows 



— I dn 

2 



2 



G a p{r,r) 

Y / dng^r) (2.24) 



In Eq. (2.24) and everywhere below the symbol J dn 
implies integration over the unit d-dimensional sphere 
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normalized to its surface area Sd 
we have e.g. 



dn = 1 



With this definition 



(2.25) 



As in the t heory of superconductivity, the solution for 
the Eq.( 2.21 ) satisfies the condition g^(r) = 1. 

Eq. ( |2.21| ) is written for a non-averaged potential u (r) . 
It is valid also in the absence of the long range potential. 
If the system is finite, Eq. (2.21) can still be obtained 



provided averaging over the spectrum is performed. 

In principle, one could average over the potential u (r) 
by expanding in this potential and averaging terms of 
the perturbation theory. This would not be considerably 
more convenient as compared to the conventional pertur- 
bation theory. Fortunately, exact solutions of Eq. ( 2.21 ) 
can be written explicitly in terms of a functional integral 
over supermatrices Q n (r) with the same structure as the 
supermatrices g n (r). 

In the next chapter we write this solution and average 
over the long range potential u (r) . 



III. SOLUTION OF THE QUASICLASSICAL 
EQUATION. THE BALLISTIC a-MODEL. 



The quasiclassical equatioBr-iEq. ( 2.21 ), has been writ- 
ten previously by MK, Ref.E2l (the long range random 
potential u (r) and the source term were not included) . 
However, they did not try to solve this equation but no- 
ticed in analogy with an equation of motion for a fer- 
romagnet that Eq. (2.21) is an extremum of a func- 
tional containing a Wess-Zumino-Novikov-Witten term. 
Assuming that this had to be a deep minimum MK 
have written the solution of the quasiclassical equation 
in terms of a functional integral with this functional. 
Unfortunately, conditions providing such a description 
have not been found as yet and the MK theory is usu- 
ally spoken of as a variational or phenomenological one. 
Although this is generally so for "equations of motion" 
like the quasiclassical Eq. (2.21), there is an important 
exception when the functional integral becomes an exact 
solution of the problem. This exception is just the case 
considered here when the equation is linear in g n (r) and 
the functions g n (r) are supermatrices. Surprisingly, this 
rather simple fact has not been reali zed b efore. 

It is clear that a solution for Eq. ( [2.21 ) is not unique 
because the equation is homogenous. To determine the 
solution uniquely one should put J ~ 0, u (r) = and 

calculate the corresponding g„ ' (r) directly without us- 
ing the quasiclassical equation for the Green function. It 
is easy to see that in this case 



ffi 0) (r) = A 



(3.1) 



Eq. (3.1) plays the role of a boundary condition. Let 
us show that the exact solution for Eq. (2.21) satisfying 
the boundary condition, Eq. (3.1) can be written as 



n (r) = Z^[J] J ^Q n (r)exp(~$j[Q n (T)]) DQ ni 



*j[Qn(r)] = Str /drdn[Af n (r)( Wi ,nV r 



-PF 1 V r u(r)V I1 )T n (r) + ( l ^ + l6 h - j( r ) ) Q n ( r) ], (3.2) 



Qn(r) = T n (r)AT n (r), T n (r)T n (r) - 1 

In Eq. ([3~2]), the partition function Z 2 [J (r)] is 



z 2 [J] 



exp 



*j[Qn(r)] )DQ a (3.3) 



and the integration is performed over the self-conjugate 
supermatrices Q n = Q n (r) satisfying the following rela- 
tions 

Ql (r) = 1 (3.4) 
everywhere in the bulk and, in addition, 

Qn± i r )\surface = Q-a± ( V )\surface (3-5) 

at the surface of the sample. 

I n or der to prove the statement that the inte gral, 
Eq.(3.2), is equal to the solution g n (r ) for Eq.( 2.20 ), we 
notice that the integration in Eq. (3.3) is performed over 



all supermatrices Q n (r) with the constraints, Eqs.(3.4 



3.5), and therefore cannot change under the following re- 
placement of the variable of the integration 



Q»(r) -» Q„(r) = U a (r)Q a (r)U a (r) 



U a (r)U a (r) = 1 
(3.6) 



On the other hand, one can formally consider the in- 



tegral, Eq.(3.3), with the transformed matrix Q n (r) as 
a functional of the transformation matrix U n (r). The 
fact that the integral Z 2 [J (r)], Eq.(3.3), does not change 
under the transformation means, in particular, that the 
first variation of the considered functional must be zero 
for any unitary matrix U a (r). 

In order to find the first variation, we make a small 
rotation, resulting in the replacement U n (r) — > U n (r) + 
<5?7 n (r){7 n (r) and compute in the linear approximation 
in the matrix (5f/ n (r) the difference between the func- 
tionals with the changed and initial matrices U n (r). 
For the supermatrix Q n (r), this means the replacement 
Qn(r) — > Qn(r) + [8Un(r), Q n (r)], which follows from 
the relation SU n (r) = — SU a (r). Then, the first varia- 
tion 5Z2 of the partition function Z%\J (r)], Eq^^) takes 
the form : 

SZ 2 [J] = -^r \ DQ n 5$[Qn] exp f-^$[Q n ] 

1 J Q2=l V 1 



(3.7) 



G 



The va riation 5<I>[Q n ] can easily be calculated from 
Eq.( |3.2|) and we write it as 

6^[Q n ] = -Str J drdn5U n (r){(v F nV 



-p^VrU (r) d n )Q n (r) + [(i^^-A - J (r)), Q n (r)]} 



+v F Str I dn (ndS) «/ n (r) Q n (r) (3.8) 
Is 



is per- 



The integration in the last term in Eq. 
formed over the surface of the sample. At the surface, 
not only Q n (r) but also SU n (r) must be invariant un- 
der the replacement n_L-^ nj^ (see Eq. ([D])) It follo ws 
from this property that the surface te rm in Eq. (13.81) is 
equal to zero. Then, substituting Eq. (3J3) into Eq. 




and taking into account that 8Z2 [J] m ust b e zero for any 
SU n (r ) we come immediately to Eq. ( |2.21 ). 

Eq.( [2.21 ) is a differential one and a class of diff eren t 
solutions ma y exis t. The fact that the integral, Eq. (3.2), 
satisfies Eq. fl2.2l|) , does not guarantee that it is equal to 
<7„ (r) , Eq. ( |2.20| ), and this should be checked separately. 
It is the supersymmetric structure of th e supermatrices 
Q n (r) that allows the integral, Eq. ( |3.2| ) , to satisfy both 
the equations. 

Putting J = an d u (r ) = one can calcul ate g n °^ (r) 
directly from Eq. ( |2.20| ) and come to Eq. (|t]). The 
s am e can be done using the functiona l int egral in Eq. 
(HH). As the functional 4>j[Q n (r)], Eq.(|3~2"|), contains at 
J = 0, u (r) = only the matrices Q n (r) and A, aver- 
aging, with such a functional gives according to general 
rulestj 



(Qn (r)>* = A 



(3.9) 




9|) with Eq. (|]l]) we conclud e tha t 
), is compatible with Eq. (2.20). 



Comparing Eqs. 
the solution, Eq 
Although Eq. (3.9) is trivially fulfilled for the superma- 
trices, it would not be necessarily correct if the matrices 
Q did not have the supersymmetric structure. In the lat- 
ter case one could obtain e.g. a non-trivial function of 
oj in the r.h.s. of it. This would invalidate th e pr esent 
approach for such symmetries. Actually, Eq. ( |3.9D tells 
us that the average density of states is a constant, which 
is the case for weakly disordered sys tem s. 

The partition functions Zi, Eqs. |J, |2.1l| , |2.12| ), 
and Z2, Eq. ( |3.3| ) are equal to unity at J = 0. As their 
logarithmic derivatives coincide for all J, we come to the 
conclusion that 



Z 1 [J]=Z 2 [J] 



(3.10) 



Thus, we rep l aced the i ntegration over electron modes, 
Eqs. (fj, |j, |2.11| , |2.12| ), by the integration over low ly- 
ing excitations that were called diffusion modes in the dif- 
fusive limit. Apparently, the name kinetic modes would 



be a more prop er on e for the limit under consideration. 
We checked Eq.( 3.10| ) additionally by a direct expa nsion 
of the both sides in J up to terms J 2 . Eq. ( |3.10] ) was 
written in Rcf.EHl at u (r) = and a — but its accuracy 
remained unclear. 

The form of the partition function Zi\J\, Eqs. (3.2, 



3.3) allows us to average it immediately over both the 
long range potential u (r) and the supermatrix M (r) . 
Assuming that fluctuations of the random potential u (r) 
are gaussian with the correlation 



(u (r)u(r')) = W(r-r'), 



(3.11) 



where W (r) is a function decaying at distances b much 
exceeding the wavelength Xp, b Xf, we average the 
partition fun ction Z <x (J) over u (r) and M (r) with the 
help of Eqs. ( 2.13 ), ( 3.11 ) and find for the final averaged 
partition function Z (a) 

Z(a)= f exp(-^[Q n (r)])£>Q n (r) (3.12) 

The "free energy" functional T takes the form 

T [Q n (r)] = T km + T imp + F^ p , (3.13) 



•TWQn(r)] = — Str J drdn[2v F AT n (r)nVT n (r) 
+2*(^^A-a)Q n (r)] 



•F imp [Qn\ 



1 I TtV 

PF 



J drdndr'dn'V l r V^(r - r') 



xStr[Af n (r) V n T n (r)]^r[AT n ,(r')V n ,T n ,(r')] 



[Qn (r) 



—^rStr ( / Q n (r)dn 



and Q n (r) = T n (r) AT n (r). 

As we see from Eqs. ( |3.12| -3.13), the free energy func- 
tional T [Q n (r)] consists of three parts. 

The first part Thin [Qn (r)] describes the kinetic modes 
in the absence of any impurities. Correlation functions of 
interest can b e obtaine d differentiating Z(d) in a as writ- 
ten in Eqs. (2.7, 2.10). The first term in Tkin [Qn (r)] 
can be written also in the form of the Wess-Zumino- 
Novikov-Witten integral with an additional variable of 
integrationt3. The second part T [Q n (r)] is responsible 
for scattering on the long range potential, whereas the 
third term is due to scattering on the short range impu- 
rities. 

Eqs. ( 3.12 - 3.13 ) are valid even in the absence of any 



impurities as long as one may use the quasiclassical ap- 
proximation. This description fails near boundaries of 
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the sample because "turning points" where the quasi- 
classical description fails are inevitable in those regions. 
Apparently, this can lead to an additional term in the free 
energy [functional analogous to the regularizer intro duced 
in Ref.lUl. In all other parts of the sample, Eqs. (3.12- 



3.13 ) are valid at all distances exceeding the wavelength 
X F . 

It is important|.tO|-epiphasize that the problem of the 
"mode locking"lljll3'Efl does not exist in the present ap- 
proach. As n 2 = 1, there are no fluctuations transverse 
to the constant energy shell and no additional averaging 
is necessary. Moreover, there are no fluctuations of the 
eigenvalues of the supermatrix Q n (r) because by con- 
struction Q n = 1. The limit of the vanishing short range 
potential r s — ► oo can be taken if some long range po- 
tential is present. The purely ballistic case without any 
randomness is more sophisticated because one needs a 
regularizer. However, this is a more delicate effect than 
the problem of the mode locking. 



The term J-i mp [Q n (r)], Eq. (3.12), does not have a 
form of a collision integral of the Boltzmann equation. 
It differs from what one obtains using the rSaddle-point 
approximation and expanding in gradients" ^l 15 !. If one 
expands the functional T [Q n (r)] in small deviations Q 
from A using e.g. a parametrization like 



T = iP + (1 - P 2 ) 



2\l/2 



PA + AP = 



(3.14) 



the first non- vanishing contribution is of the order of P 4 , 
which means that Ti mv [Q n (r)] does not contribute to 
the bare propagator at all. .— . 

At the same time, it was demonstrated!!!!! that the bal- 
listic CT-model with a term in the form of the collision 
integral could be obtained at large distances exceeding 
a single particle mean free path I. In the next chapter 
we will clarify this question by demonstrating that the 
standard form of the collision integral can really be ob- 
tained at large distances as a result of a course-graining 
procedure. 



IV. REDUCED BALLISTIC a-MODEL 



The free energy functional T [Q n (r)], Eqs. ( |33^ - pl3| ) 
is most general and applicable at all distances exceeding 
the wavelength Xp- We neglect now the part T s [Qn (r)] 
originating from the short range impurities and concen- 
trate on studying properties of the long range scattering. 

At very large distances exceeding th e mea n free path l tr 
(we evaluate this length later, see Eq.( 4.3Qb ) the er-modcl 
must acquire the standard diffusive forrrB. However, if 
the long range potential u (r) is weak such that l tr 3> fe, 
one more intermediate scale is important for describing 
the behavior of the system, namely, the one related to the 
Lapunov exponent Xl = T i}- The corresponding time tl 
is a time required for two particles moving initially par- 
allel to each other to increase the distance between them 
by a factor of order unity (This definition gives, of course, 



the order of magnitude of tl only). The importance of 
this time was first pointed out in Ref.EZl, where a weak 
localization correction was calculated for scattering on a 
long range potential. This time is also releptnt for corre- 
lations of wave functions discussed recently^ for a model 
of weak scatterers where it was demonstrated that only 
at distances R exceeding the II = vftl the notion of 
separate diffusons makes a sense. 

The length II was estimated for the model of long 
range weak scatterers corresponding to the case consid- 
ered here as 



h 



hr (b/krf 3 



(4.1) 



which shows that II is between b and lt r , b <C II <S hr- 

At distances R smaller than II, two particles that 
started their motion along parallel trajectories still move 
parallel to each other and, following Ref.tZl, we call the 
region inside II Lapunov region. 

In this chapter, we wan t to show that the ballistic 
non-linear g-model, Eqs. (3.12, 3.13), acquires a more 
familiar formaa provided an integration over Q n (r) 
within the Lapunov region is performed. This integra- 
tion is some kind of the course-graining procedure used 
very often in statistical physics. As a result of such an 
integration, we obtain another field theory that can be 
called a reduced ballistic er-modcl. 

Proceeding in the standard way we separate fluctua- 
tions of the matrix Q n (r) into fast and slow parts. By 
fast variations of the supermatrix Q n (r) we mean fluctu- 
ations changing fast within the Lapunov region and the 
rest is classified as slow ones. This procedure is very sim- 
ilar to the renormalization group (RG) scheme applied in 
2D in the diffusive region (see, e.gtl). Carrying out the 
course-graining procedure when deriving the RG equa- 
tions is simplified by the fact that all arising integrals 
are logarithmic and one needs to know only the order of 
magnitude of the cutoffs at each step of the shell inte- 
gration. In the course-graining procedure used here, re- 
sulting integrals are not logarithmic. At first glance, this 
would make the entire procedure rather tricky because a 
renormalizcd free energy functional would depend on the 
length I l estimated by the order of magnitude only. 

Fortunately, there exists a method of integration over 
the fast modes resul ting in a renormalized free energy 
functional P[Q n ] Eq.(4.22) that does not contain II as a 
parameter. This length will be implied for P[Q n ] as an 
ulraviolet cutoff only. 

The separation into slow and fast fluctuating parts is 
performed as follows 



TJj) = T n (r)V n (r) 



(4.2) 



where both P n (r) are V^(r) are unitary supcrmatrices. 
The supermatrix T„(r) is supposed to describe slow 
modes, whereas Vn(r) is res ponsi ble for fast ones. Substi- 
tution of Eq.(fO) into Eq.(3.13) results in a new action 
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in which the two fields interact in complicated way. We 
write this action as follows 



T km [T}=T km [V}+T' km [f,V] 
J~imp [T] = J~i, np [V] + J- imp [T, V] 



where 

T kin [V] 



— Str I drdn[2v F AV n (r)nVV n (r) 

- 2mA - 2\ l )AQW ( r )] (4.4) 



— Str I drdn[2v F Q^f n (r) nVT n (r) 

+i(f n (r) (u + iS - 2aA)AT n (r) 

-(u + i5)A + 2a)Q^(r)] (4.5) 



:5tr[AV n (r)Viyn(r)]S , tr[AV I1 ,(r')Vi^ n /(r')] (4.6) 



nnJr, V] = ~ f — ) / drdr'dndn'Vj.V^ (r - r') 



PfJ 



x [^(Q^ (r) & n (v))Str(Q^ (r') (r')) 



+2Str(AV n (r) ^MStr^) (r) <F n , (r'))] (4.7) 



In Eqs. (4J-4/7 ) 



g(°»=l/ n (r)Ay n (r) : 



$1 (r)=f n (r)V 4 „f n (r) 



(4.8) 



The separation into the fast and slow parts, Eq. (4.2), 
is not simple for the problem involved. This is because, 
for a given n, a characteristic dependence of Q n (r) on the 
coordinate r is extremely anisotropic and one cannot in- 
troduce an isotropic momentum shell for integration over 
the fast modes. Besides, we should not violate the rota- 
tional invariance when integrating over the fast modes. 
This goal is achieved by writing th e ter m with the La- 
punov exponent Xl in .Ffcm[V], Eq. (4.4). This is similar 
to an invariant integration over a momentum shell for a 
diffusive <7-model. 

Introducing a notation 

(...}o = J (...)exp(-F kin [V])DQW 
we write the reduced free energy F as 



F = -In(exp(-^ in [f,Vl- 



(4.9) 



(4-3) (due to the supersymmetry (1)q = 1] 



We calculate the average in Eq.(4.9) by expansion of 
the exponential and computing averages of all terms ob- 
tained in this way. In the first order, using the relation 



(QL 0) (r))o = A 

we obtain 

(T' kin [T,V]) Q =T km [Q n (r)] 
with T km [Q n (r)] from Eq. ( fig ), 

{Fimp [V]) = 0, 



(4.10) 
(4.11) 
(4.12) 



which is due to the supe rsym metry, as well as the second 
term in ^„„[f , V], Eq.Q, 



WL.Ef.vi).— i(= 



drdr'dndn'Vi V J T ,W (r - r') 



x (Str(Q^ (r) < (r))Str(Q^ (r') (r')))o (4.13) 

If we re placed the average of the product of Qn°' in 
Eq. ( 4.13 ) by the product of the averages we would sim- 
ply come back using Eq. ( 4.1C ) to Ti mp [Q n (r)] in Eq. 
( |3.13| ). However, now the supermatrices Q n (r) vary at 
much longer distances than the radius b of the corre lation 
function W (r — r'). As the integrand in Eq. ( |4.13| ) con- 
tains derivatives V r of W (r — r'), the integral is small 
and can be neglected. 

The main contribution comes from the irreducible part 
of the correlation function which we denote as 

(((g(°)(r)) Q/3 (g^(r')) W ))o (4.14) 

(at the moment we do not distinguish between the ad- 
vanced/retarded blocks and the others and write all in- 
dices standing for matrix elements as superscripts). We 
calculate this average in Chapter V in a gene ral case in- 
cluding the impurity potential (see Eq. |5.9| ). Now we 
neglect the impurity potential and make the replacement 
6 — > X l when determining the fu nctio n Q nrl > (p) from 
Eq.(|J). Then, we substitute Eq.Q into Eq.(|3|) 
and transform the result of the substitution by introduc- 
ing new integration variables R = , p = r — r'. Then, 
we note that the integrand contains the function Q nn ' (p) 
that changes on the distance II (the limit lu <C Xl is 
implied), while the supermatrix$^ (r) changes slower. 
This allows us to rewrite Eq. ( 4.1 3| ) as follows 



2pl 



Str / drdndn' 



x«'( r )) X / d P V i W j W(p)Ag an ,(p) 



(4.15) 
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The notation (• • ■) means here the part of the superma- 
trices anticommuting with A. Eq. ( 4.15 ) can be reduced 
to a more simple form. In the Fourier trans forme d rep- 
resentation, the integral 7 nn ' over p in Eq. ( 1.15 ) takes 
the form 



/ ^ W(q)<w(-q) 



(4.16) 



Solving Eq.(5.8) we find for the Green function in the 
limit lu <C \l 



£nn'(-q) 



WFnq + AlA 



(4.17) 



The function W (q) decays at momenta much smaller 
than pp and only such momenta give the main contri- 
bution in Eq. ( 4.1 6] ) . Therefore, we can appro xima tely 
write the function Gnn'(— q) in the integral, Eq.(4.16), as 



iSr, 



(2*) 



2m 



(4.18) 



where p = pi?n. 
Using the fact that 



we reduce Eq. ( |4.15| ) to the form 



^L P [T, V]j Q = ^Str J dndrV^Q (r) V J n Q (r) (4.19) 



(27T)' 



2m 



Cha nging the variables of the integration in the integral 
Eq.( 4.19 ) to p' = p q and integrating separately over 
£ = p' 2 /2m — £f and n' = p'/p' we obtain finally 



T'- \f V 



dr / dndn'(n — n')'(n — n 



(to/) 



:W nn / Str^Qn (r)ViQ n (r)] (4.20) 



where W nn > = W (pf (n — n')). This result corresponds 
to the collision integral in the Boltzmann kinetic equation 
in the limit of small angles of the scattering. 

Proceeding further we can calculate, in principle, not 
only (T' imp \T, V])q but also its higher cummulants like 



To estimate these cummulants we should consider an av- 
erage of the type 



((VnVn'QW (r)Q<°> (r') 
xV ni V ni g( 1 ° 1 ) (r 1 )Q[ i °, ) (r' 1 )»o 



(4.21) 



where r — r'| ~ |ri — r[\ < b and we do not write explic- 
itly indices. 

Th e cor relation function for the product of four Q^°\ 
Eq. ( 4-21 ), resembles a correlation function calculated 
in RefJi3 where it was demonstrated that, being a com- 
plicated function of small distances, the reducible cor- 
relation function decoupled into 2 diffusons as soon as 
the distance between the points r, r' and ri,r[ exceeded 
the Lapunov length If- This means that the irreducible 
correlation function of the type as in Eq. (4.21) decays 
at distances of the order of the length If- As the sec- 
ond cummula nt c ontains a product of four slow functions 
$ n (r), Eq. (18), we conclude that the cummulant ex- 
pansion is effectively a series in ^V r . Therefore, beyond 
the Lapunov region we may keep the average ^ mp [T, V] 
only. Using Eqs. (4.9- 4.12 , 4.20 ) we write the reduced 
ballistic cr-model applicable beyond the Lapunov region 
in the form 



F[Q n (r)] = —Str I ,/r 



{2v F AT n (r) V r T n (r) 
) n (r))dn 
lU nn , (n - n') 1 (n - n') J 



+£(u> +iS- 2aA)AQ n (r))dn 



xVj 1 Q n (r)V J n Q n (r)dndn' 



(4.22) 



The reduced ballistic cr-model, Eq. ( 4.22j ), has the same 
form as the cepcesponding cr-model derived for a long 
range potentials provided the limit of small angle scat- 
tering is considered. The non-linear cr-model was ob- 
tained in Ref.E3 under the assumption that all distances 
exceeded the single particle mean free path I using, as the 
first step, the saddle point approximation. However, the 
saddle point approximation is not good for a long range 
potential and, hence, the length / cannot be a good quan- 
tity. Now we see that one should use the length If- The 
reduced cr-model is applicable at distances larger than 
If and this agrees with the suggestion of Refill. It is 
remarkable, that the reduced cr-model does not contain 
explicitly the length If as a parameter. It can enter as an 
ultraviolet cutoff only. Therefore, it is sufficient to know 
If by order of magnitude. 

On e can further simplify the reduced ballistic cr-model, 
Eq. (4.22), at distances exceeding the transpor-Lm^aa 
free path l tr 3> If- This route is well developedEJ'EJEj. 
Separating again fast modes (strongly varying at dis- 
tances smaller than l tr ) from slow ones with the zero an- 
gular harmonics (fluctuating at distanc es exc eeding l tr ) 
we write the supermatrix T n (r) in Eq. ( 4.22 ) as 



T n (r) = U(r)TW ( r ) 



Then, the free energy functional F[Q n (r) 
takes the form 



(4.23) 



Eq. (4.22) 
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F[U,T<®] = F [tW]+F'[U,tW], 



(4.24) 



Fo[T (o)] = ™str J [2v F T^ (r)nV r T(°) (r) 

+i(u + iS)AQ^(r)]de (4.25) 



F' 



UjT (0) 



Str I dr[2v F Q ( n } (r) n$ (r) (4.26) 



717/ 

T 



+^ (w + i<5) (Z7 (r) At/ (r) - A) Q^ 0) (r) + 
| W nn , (n - n') 1 (n - n') J (r)V> B QW (r)dndn'] 



where QL 0) (r) = T^ 0) (r) AT^ 0) (r), $ (r) = 17 (r) VU (r). 

The next step is to integrate over the fast modes, which 
leads to the free energy functional F<nf [Q] 



F dif [Q] = (exp(-F' 



n T (o) 



(4.27) 



where 



<-. >o= / (...)exp(-F [rW])DQ 



)(0) 



Integration over (r) can be performed using, e.g. 

the parametrization 



(r)=A(l + iP n (r))(l-iP n (r))- 



(4.28) 



with P n = -P n , P n A + AP n = 0. 

Then, one can expand P n (r) (as a functio n of n ) in 
angular harmonics an d exp and in P n in Eqs. (4.25). In 
the first term in Eq. fl4.26| ), it is sufficient to keep only 
the linear term in P n , while in the third term we should 
expand up to quadratic terms. Due to the presence of the 
linear term integration over the first harmonics is most 
important. Performing this gaussian integration we come 
to the diffusive cr-model 

Fdif [Q] = Y str fl D ( V( 5) 2 + 2l ^ + A ® dr 

(4.29) 

where Q (r) = U (r) AU (r) and D = v F r tr /d, d is the 
dimensionality. 

The transport time T tr can be written through the 
function W nn > as 



I 1 = 2tiv J W nn < (1 - nn') da! 



(4.30) 



whereas the transport mean free path l tr is h r = i>FTt r - 
If in addition the short range potential is present, the 
effective mean free time r e // can be written as T~fj = 



Writing the reduced cr-model, Eq. ( 1.22] ), in a limited 
volume one can come rather easily to the Wigner-Dyson 
level- level correlation functionsB. This can be done sep- 
arating fluctuations of T n with the zero space and angle 
harmonics from other degrees of free dom . We write this 

" HI) 



separation in a form similar to Eqs. ( 4.2 



T n (r) = UT^(r) 



(4.31) 



where U depends neither on the vector n nor on the coor- 
dinate r, while T„ (r) contains only non-zero harmonics in 
the both coordinates. Integration over T„ (r) can b e per - 
formed using a parametrization analogous to Eq. (4.28). 
Excitations corresponding to the matrix T' n (r) have a 
gap and their contribution can be neglected in the colli- 
sion region in the main approximation in the parameter 
cumin far , tl) ■ Then, we come to the zero-dimensional 
version F [Q] of the c-modcl 



Fo[Q] = 



7ri (to + iS) 
4A 



Str (AQ) 



(4.32) 



(vV) (V is the volume) is the mean 
This leads directly to the Wigner-Dyson 



where A 
level spacing, 
statistics! 

The above discussion describes completely the method 
of derivation of the cr-model based on quasiclassical 
Green functions in all regions. However, explicit calcula- 
tions with the ballistic cr-model are not simple and cal- 
culational schemes are not necessarily the same as those 
used in the diffusive limit. The most unusual is the La- 
punov region. In the next chapters we consider calcula- 
tional schemes that can be useful for clean systems. 



V. CORRELATION FUNCTIONS FOR A LONG 
RANGE DISORDER 

The ballistic cr-model derived in the previous chap- 
ters has a form of a functional integral over the super- 
matrix Q n (r) with the constraint Q\ (r) = 1 . Using 



parametrizations of the type of Eqs. (3.14, 4.28) and ex- 
panding in P n is a potentially dangerous procedure be- 
cause the rotationalrinvariance may be lost. From study 
of disordered metalsH we know that such a perturbation 
theory works well in the diffusive region but fails in the 
regime of strong localization. Interestingly enough, the 
perturbation theory in terms of the ballistic excitations 
does not work also in the Lapunov region. Using the 
expansion in P n one comes to propagators 



(u^nq + loA)~ 



(5.1) 



and integrals over both n and vectors q. 

These integrals are not generally convergent and the. 
result depends strongly on a regularization. In Ref.EHl 
the regularization was carried out by writing an addi- 
tional term in the cr-model that would arise if a small 
amount of short range impurities was added. It was 
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suggested that, at the end of calculations, the regular- 
izer ,-Gpuld be put to zero. In contrast, the authors of 
Ref.tZl assumed that a regularizer had to remain finite 
and demonstrated in a subsequent publication^ that its 
presence in the cr-model leads to an anomalous contribu- 
tion. Problems with diagrammatic expansions, in the bal- 
listic excitations were also discussed in RefEd. The fact 
that one can hardly speak about separate diffusons (kiiie* 
tons) in the Lapunov region has been noticed in Ref.tZl 
and emphasized in Ref.Eil 

At the same time, the integration over Q n in the ballis- 
tic cr-model, Eq. ( |3.12| , 3.13) is well defined and the prob- 
lems aris e only afte r using a parametrization like those 
in Eqs. ( 3.14 , 1.2S| ) and a subsequent expansion in P n . 
Therefore, we should understand what one can do if one 
may not use such a parametrization. It turns out that 
the most proper way of computations of correlation func- 
tions is deriving equations for them. This is analogous to 
what one does in models for turbulenceE3. Unfortunately, 
this method is not as general as the perturbation theory 
but for some correlation functions closed equations can 
be derived without difficulties. 

As an example, we consider the function y 00 (ri, r 2 ; u>), 
Eq. (|]^) (without averaging over the long range poten- 
tial). Using Eqs. ( pT0| , pT0| ) we reduce this correlation 
function to the form 



y°°(n,r 



2< 1 



2(7T^) / dnidri2 



^ 2 _ i Q 8 n 4 1 (ri)QS(r 2 )exp(-^$ [Q n(r )])^Q n 

(5.2) 



with the functional <f>o[Qn(r)] from Eq. (3.2) in which 
we neglect the short range impurities and put J = 0. 

In order to find the average (<5„t ( r i)Qn2 ( r 2)}<2 (here 
and below the symbol (. ■ -)q is used for averaging with 
the functional < &o[ ( 9n( r )]) we introduce the function 
5n (r; a n ) (superscripts are omitted): 

g n {r;a n ) = Z^i&n] / Q n (r)exp (-^-<f>a n [Qn}) DQ r 
JqI=i v 1 ' 

(5.3) 

Z [a n ] = [ exp (-^aJQn}) DQ n 
JQI = 1 X 2 ' 

$a„[Qn] - $o[Qn] - iStr J drdna n (r)Q n (r) (5.4) 
Below, the source d n is assumed to be a matrix anti- 



longer self-conjugate: a„(r) ^ a n (r). Then, the first or- 
der c/n 1 ^ (r; a n ) of the expansion of the function g n (r; a n ) 
in the source gives the irreduceable correlation function 



commuting with A. Besides, in contrast to Eq.(3.2), it 
includes now the dependence on the direction n and is no 



(r;a n ) =i^-J dr'dn' «Q n (r) 

xStr(Q n , (v')a n , (r'))»Q 



(5.5) 



At the same time, the function g n (r; q n ), Eq .(|5.4|), is a 
solution for the equation (c.f. with Eq. ( |2.21 )) 



(u_FnVr - Pp 1 V r u(r)3„) g n (r; a n ) + 

i(lu + iS) 

'- [A, 5n (r;a n )] = 

7 _ 

- [a„ (r) + a n (r) , g n (r; a n )] 



(5.6) 



Eq. fl5.6|) can be solved by iterations expanding in a n (r). 
As the zero order in a n (r) we may take g^ (r) = A using 
again the supersymmetry of the integral Eq.( |5.3| ) when 
the source is disregarded. Then, we obtain in the first 
order 

<£> ( r ;a n ) = 1 J dr'dn'(g nD , (r,r') 

+kg nn , (r, r') A) [a n , (r') + ~a n (r) , A] (5.7) 

The kernel £?nn'(r,r') in this expression is a Green func- 
tion for the equation: 



+i(w + i(5)Ag nn /(r,r') = <WJ(r-r') 
iFrom Eqs.(|J-|jl) we find easily 



(5.8) 



1 

2m> 



<C(r,r')A^+ 



(((Qn(r)) /3Q (Q n -(r')) W )) Q 

+ A0 "GZ> ( r , r') A'" - (r, r') A^<F Q - 
-A Pv GZ> (r, r')<F a ] + c.c. (5.9) 

where c.c. st ands for "charge conjugated" terms. Then, 
taking in Eq.(5.9) different values of the superscripts one 
can find corresponding integrals. For example, putting 
a = 7 = 4, (3 = n = 8 results in the relation 



and we obtain finally 

r 00 (ri,r 2 ,w) = 4-7TI/ J dn 1 dn 2 V(l;2) 



(5.10) 



(5.11) 



The function V(l; 2) used in Eqs. ([5.10|]5Tl| ) is the Green 
function of the Liouville operator C — — -Jr- and 

^ op Or Or Op 

satisfies the equation 

(w F niV ri - p^VrM^d^ - i (w + iS)) V (1; 2) = S h2 

(5.12) 
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We used here the notations j = (nj,rj), j = (— nj,rj) of 
RefB. 

The quantity Y~ 00 (ri, r 2 , to), Eq.(| is an extension 
of the non-averaged density-density correlation function 
to arbitrary scales. We emphasize that the result is exact 
within the quasiclassical approximation. 

A similar computation can be carried out for a higher 
order correlation function Y 2 00 (ri, r 2 , f3, r 4 ; us) 



Y 2 00 (ri,ra,r3,r 4 ;w) 



(5.13) 



= Gf (n, r 2 ) Gt„ (ra, r x ) G? (r 3 ,r 4 ) G^_ w (r 4 ,r 3 ) (5.14) 

This function can be written as an integral over the su- 
pervectors ip (r) as 



F 2 00 (n^.r^jw) 
= 16 f V 4 (r 2 )^ 4 (ri)V 8 (ri)^ 8 (r 2 ) 
x i> 2 (r 4 )^ 2 (r 3 )^ 6 (r 3 )^ 6 (r 4 ) e - L W 



(5.15) 



(5.16) 



Using the res ults o f Chap. Ill we rewrite the integral 
over i/j (r), Eq. ( |5.15| ), in terms of an integral over Q n (r) 
as 

4 r 

r 2 00 (r 1 ,r 2 ,r3,r 4 ;w) = (2tt I ,) 2 [] / dmM(l, 4; 2, 3) (5.17) 

i= 1 

7W(1,4;2,3) = 



(?) 



^ ) «£(l)«£MQ5(*»)Q£( r «)>9 (5.18) 



I n or der to calculate the corr elat io n fu nction Ai, Eq. 
( 5.18 ), one should expand Eqs.(|5.3|), (5.6) up to the third 
order in the source and then compare them with each 
other. Writing the integral Eq.(5.3) for a — 8, (3 = 4 and 
assuming that only the elements a^(r), a^ 6 (r), d„ 2 (r) 
of the source are not equal to zero we expand it in the 
element a^ 8 (r). Keeping only the first order in the ex- 
pansion and putting a„ 8 (r) = ti5 nin ,i5( ri — r 2 ) we obtain 
the following integral generalizing Eq.(5.10) 



X> nin2 (ri,r 2 ;d^) 



~2 



xe 



-*a' [Q»W] 



DQ n 



(5.19) 



The new source a' n in Eq. ( 5.19| ) differs from the previ- 



ous one a„ by the substitution a^ 2 (Y2) = 0. At the same 
time, making the same for the solution of Eq.(| 
obtain 



we 



[w F n 2 V 2 - p F 1 \7 r2 u(r 2 )d n2 
+i{uj + iS)]V nin2 (n, r 2 ; a' n ) 



1 



<5i2[5 88 2 (r 2 ;0-^4( r2 . 0]Z2K] 



(5.20) 



In order to come now to the integral, Eq.(5.18), one 
should find the term of the expansion V niIl2 (ri, r 2 ; a' n ) in 
the source a' n bilinear in the elements d„ 6 (r), d^ 2 (r). As 



concerns the elements g„ 2 (r 2 ;dn), g 8 1 8 (r 2 ; a' n ) , they sat- 
isfy Eq. ( |5.6|) implying the replacement d n (r) by a' a (r). 
Due to the structure of the source each of the equations is 
closed and contains no terms with the source. Their so- 
lutions are therefore ±1, respectively. Expanding further 
Z2[a' n ] in the elem ents d„ 6 (r), a„ 2 (r) and substituting the 
result in Eq.( 5.20| ) we obtain 



[u F n 2 V 2 - p F V r2 w(r 2 )5 n2 

-t(w + iS)]M(l, 4; 2, 3) = <5 12 P(3; 4) 



(5.21) 



The operator in the l.h.s. of Eq.( [5.21 ) acts on the vari- 
ables of a one of the "diffusons" . In order to get also an 
equation with respect to the varia bles of the other diffu- 
son we consider the integral, Eq.( [5.3|) , with the source 
having now only the elements a£ 2 (r) , d^ 8 (r), d 84 (r). 
Putting there a — 6, (3 — 2 and repea ti ng su ccessively 
all the steps that lead us to Eqs.(5.20)-( |5.21 ) we come 
to an equation like Eq.( 5.21 ) in which n 2 , r 2 in the left 
side are replaced by n 3 , r 3 and the variables 1, 2 in the 
right side by 3, 4 respectively. Adding this equation with 
Eq.(|5 21D we find 



(w F [n 2 V r2 +n 3 V r3 ] -p F 1 [V r2 M(r 2 )a„ 2 + V r3 ii(r 3 )<9 n3 ] 
-2i(w + iS))M(l, 4; 2, 3) = 5^(3; 4) + 8 34 T>(1; 2) (5.22) 

Eft. ( [3.22] ) agrees with the corresponding equation of 
Ref.U. This equation is written before averaging over the 
long range potential and is again exact within the qua- 
siclassical approximation. However, the averaging over 
the long range potential is_uot trivial. Analyzing this 
equation Aleiner and LarkirO demonstrated that the av- 
eraged function 7W(1, 4; 2, 3) decouples into 2 diffusons 
(one can connect the Green functions G R ' A in 2 differ- 
ent ways forming 2 diffusons) at lengths exceeding the 
Lapunov length 1^. We call here this region collisional. 
A similar analysis starting from a different formulation 
was carried our in Rcf.E3 where the authors came to the 
same conclusion. This corresponds to the possibility of 
expanding in small fluctuations of the supermatrixQn 
usi ng a p arametrization like the one given by Eqs. ( 3.14 ) 
or $L2% . 

In the Lapunov region, such an expansion is impossible 
and the only way to analyze correlation functions is to 
write equations for them for a fixed potential and then 
average the solution over the potential. Unfortunately, 
we do not know how to derive exact equations for av- 
eraged correlat ion functio ns. The form of the ballistic 
a- model, Eq. ( 3.12 , |3.13| ), does not allow us to derive 
such equations because of the complicated form of the 
teim J^imp' 

Apparently, using the ballistic cr-modcl for calculations 
in the Lapunov region does not bring considerable advan- 
tages. At the same time, we see that the fi-model is ap- 
plicable also in this region, although one cannot use the 
perturbation theory. Therefore, the conclusion of Ref.Ej 
that a completely different theory should be constructed 
for the Lapunov region is too pessimistic. 
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VI. PERIODIC ORBITS 

In this chapter, we consider a clean finite system (quan- 
tum billiard). In principle, there can be holes (antidots) 
in the system but this can be discussed in terms of a more 
complicated surface. 

A standard tool for computation of the one particle 
density of statea-is, the Gutzwiller trace formulae! (for a 
review, see, e.g.tmi]). This formula allows one to express 
the non-averaged density of states in terms of a sum over 
periodic orbits. This sum is actually divergent but there 
are methods to obtain reasonable results from it. As long 
as energies involved are of the order of the inverse period 
of short orbits one can extract a detailed information and 
compute also level-level correlation functions. Very often 
a statistical information is of the main interest and one 
calculates quantities averaged over spectrum. 

The situation becomes considerably more difficult if 
one studies behavior at small energies of the order of the 
mean level spacing A. Even calculation of the first non- 
vanishing non-oscillating terms for the Wigner-Dyson 
statistics is not pimpl-jr 3 ! and difficulties grow when cal- 
culating next ordersc-i On the other hand, the Wigner- 
Dyson statistics can easily be obtained from the zero- 
dimensional cr-model and the only question is when this 
OD cr-model description is valid. 

In Chap. IV we came to the conclusion that the OD cr- 
model can be used for the model of weak long range scat- 
terers if the frequency uj is smaller than min (t^ 1 , 7 y~ ) 



For the level-level correlation function R(uj), Eq. (2.7), 
the parameter w is the energy difference between two 
levels but what are the parameters Tt r and tl or, respec- 
tively, l tr = v F T tr and l L = VfTlI 

For a clean quantum billiard the transport mean free 
path should be of order of the system size L. At the 
same time, an important role should be played by the 
Ehrenfest time 



t E = A" 1 ln(L/A F ) 



(6.1) 



where A is the Lapunov exponent for scattering on the 
boundaries. The corresponding length Ie = vptE is much 
larger or of the order of the system size L. This time de- 
termines t.Vip crossover from the classical to the quantum 
regimeEjclj. The question about the Ehnjpfest time is 
becoming popular in mesoscopic physicsE3n3nJ. In the 
present work, we are not able to obtain this time within 
the approximations used. Averaging over the spectrum 
is equivalent to averagin g over infi nite range impurities 
and we see from Eqs. ( 3.12 , 3.13 ) that only the term 
^ r fcin[Qn(r)] is present in this case. 

At first glance, it was not necessary to average over the 
energy when deriving the quasiclassical equations and the 
ballistic cr-model in Chap. II. However, this depends on 
what limit is taken first: the infinite size of the systems 
at finite disorder or vanishing disorder in a finite system. 
In the former case, an additional averaging over the en- 
ergy is really not necessary. However, energy levels of a 



finite system are quantized and we cannot directly fol- 
low the arguments of Chap. II for the latter c ase. F or 
example, the function g n (r) introduced in Eq. ( 2.2C ) is 
not a smooth function because one should sum over £ in- 
stead of integrating over it. However, the averaging over 
the energy improves the situation and makes possible us- 
ing the quasiclassical equations. Since we sho uld a verage 
the partition function with the sources, Eq. (2.11), and, 
hence, the Green functions, one can simply add aver- 
aging over the energy to the summation over £ in the 
definition of g n (r), Eq. (2.20). Then, the function g n (r) 
is a smooth function of r and we can repeat all the 
subsequent arguments leading to the cr-model. So, al- 
though there is no disorder in the free energy functional 
Fkin [Qn (*■)], the energy averaging is implied for the clean 
quantum billiards. 

As we have mentioned, our quasiclassical approach 
should be valid everywhere except in the vicinity of the 
boundaries, where the approximation fails near the turn- 
ing points. Therefore, a more accurate computation 
might produce an additional term in tho-ballistic cr-model 
near the boundary. Aleiner and LarkinO introduced such 
a term modelling the quantum diffraction by fictitious 
short range impurities. 

We think that an additional term in the cr-model due to 
the quantum diffraction is really necessary but leave its 
derivation for a future work. Instead, we will try now to 
derive the function R (uj) neglecting this term. In other 
words, we consider energies uj exceeding t^ 1 whatever it 

In RefM an additional term was added as a regular- 
izer, which had to be put to zero at the end of the cal- 
culations. However, proceeding in this way the authors 
of Ref.E3 got a result that did not agree with the one 
obtained from the Gutzwiller trace formulate. The dis- 
crepancy has been called "repetition problem" and was 
discussed in a number of worksEJiij. 

We want to show now that the result obtained with the 
ballis tic cr- model containing only the term JFfci n [Q n (r)], 
Eq. ( 3.13 ), agrees with what one can expect from the 
trace fptmulae. In contrast, the perturbative approach 
of Ref.fca is not accurate in this limit. 

We start our discussion with Eqs. ( 3.12 , [3.13| ), where 
only the term ■Ffci t JQ n (r) | is left. The level-level correla- 
tion function R(u>), Eq. fl2.7|) , can be written as 



R(u) = l- Rel(w) . 



(6.2) 



K ' 2V 2 



x {Q 8 J (r') + l)) kin drdr'dndn' (6.3) 



where we use Eqs. (6/7), ( p.8[) and introduce a notation 
< ... >kin= f (...)exp(-^[Qn(r)]) DQ n 
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T he fre e energy functional ^F^ n is obtained from T\ 



kin i 



Eq.(p^|) by putting a = 0. 

Due to the absence of any regularizer a perturbation 
expansion in ballistic excitations (diffusion modes) can- 
not be good because one obtains d iver ging integrals with 
propagators like the one in Eq. (5.1) in the integrand. 
Therefore, this method should not be applied and we 
should try something different. At the same time, the 
specific form of the functional ^Sl[Qn ( r )] that does not 
contain second space derivatives allows us to simplify the 
functional integral by reducing it to functional integrals 
on periodic orbits. 

In order to proceed in this way we discretize the phase 
space writing the functional integral as a definite integral 
over Q at all sites of a lattice in the p hase space. In this 
way we write the function I (w), Eq. (3.3), as 



n«) = [ ^§£ E <(<#('<)-i) 



x(Q?>i) + l))w 



(6.4) 



The free energy functional on this lattice takes 

the form 



4l{Qn(r)]=TF n AQ} 



(6.5) 



In order to reduce the multiple integral over all Q on 
the lattice sites in the phase space to a simpler form we 
use the following equalities that can be proven using the 
methods of integration over supermatricesEl 



exp(-F n [Q]) J DQ = l 



J (Qi 4 (r) - 1) cxp (-F n [Q])DQ 
= J (Q« 8 (r) + l)exp(-F n [Q])/JQ = 



(6.7) 



(6.8) 



for any r and n. 

Let us understand first how to simplify the function 



I (u>), Eq. (S.3), for an infinite sample. Using Eqs. (3.7 



6.8) we conclude immediately that on ly t he terms with 
n.; = iij contribute in the sum in Eq. (6^4). As concerns 
the free energy function al, o ne can integrate over all Q n 
with iij =/= iii using Eq. (6/7), which leaves only one term 
F ni [Q] in the exponential. Only the term with parallel 
to the line connecting the points r and r' gives a non-zero 
contribution and, moreover, integration over Q on sites 
outside the line gives 1 because the free energy functional 



J^[Q n (r)], Eq. fl6.5|), does not contain couplings of Q 



on these sites with Q on sites on the line between r and 



F n [Q] = \ U Str[ E -pA^n (r 4 ) T n (r{) 



-^i{uj + i8)KQ n (r;)] 



In Eqs. (3. 4, |6.5| ), / is the elementary length in the 
coordinate space along a trajectory and [AO] is the ele- 
mentary phase volume in the phase space perpendicular 
to it (i.e. to the unit vector n). The summation in the 
first term in -F n [Q] is performed over nearest neighbors 
on the trajectory and in a certain order. The elementary 
volume [AO] in Eqs. (6.4, |6.5|) can be written as 



d-l 



[AO] = S^ 1 J] (AniAr^ 



(6.6) 



where Si is the surface of th e unit sphere in the d- 



dimensional space. In Eq. (3.6), Ar l ± and Arij_ are the 



elementary length and momentum in a direction perpen- 
dicular to the path. 

In principle, the length / and the space volume [AO] 
may be arbitrary. At the same time, we should remem- 
ber that we have used the quasiclassical approximation 
and the length / may not be smaller than the wavelength 
Xp. However, only for specific choice of [AO], the func- 
tional F n [Q], Eq. (6.5), remains single valued. As we 
will see, this choice corresponds to the Bohr-Sommerfeld 
quantization rules. 



Thus, we come to an integral over all Q on sites along 
the line connecting the points r and r' 



I(u)=I+(u)+I-{u), 



(6.9) 



i± M = [ ^§f E / (Qf (*0 - 1) (Qf to) + 1) 

xexp(-F + [Q]-F4Q])DQ, 



F ±[Q±] ^^mistr[ £ ^AT ±(ri )T ± (rO 



{ r i> r S}± 



/ 



+ E*(^ + ^) A( 9± (n)] 



The signs + and — correspond to different directions of 
the trajectory; the summation over the pairs {ri,r-}-i- 
in the free energies F± [Q±] should be fulfilled in the or- 
der conforming with its direction. One can see that the 
free energies and, hence, both the terms I + (ui), /_ (u>) 
are equal to each other. This relation follows formally 
from the definition of the conjugation and the equality 
Qn( r ) = Q n ( r )- ^ is important that the time-reversal 
symmetry is not violated. 

We see from Eqs. ( |6.9[ ) that the functional integral, 
Eq. (6.3), over Q a (r) on all sites in the phase space of 
an infinite sample has been reduced to a functional inte- 
gral along a line and averaging over all directions of this 
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line. This is due to a specific form of the free energy func- 
tional containing only first space derivatives. Any regu- 
larizer containing second derivatives would make such a 
reduction impossible. 

What happens if the sample is finite ? W e can reduce 
as before the functional integral, Eq. (|6.3| ), to an inte- 
gral over the line. However, we can follow this line until 
we reach the boundary, where we have a degeneracy that 
follows from the boundary condition, Eq. (3.5). Namely, 
the supermatrix Q n belongs to 2 different lines. This 
means that, having reached the surface, we can follow 
the line obtained from the first one by a specular reflec- 
tion. We can keep going along the second line until we 
reach another boundary, etc. In principle, we have two 
possibilities: 

1. After several reflections from the boundaries we 
come to the same point in the phase space or, in other 
words we get a periodic orbit. 

2. The broken line obtained after the reflections on the 
boundaries does not close in the phase space, which can 
be considered as a periodic orbit with an infinite period. 
Of course, we should speak rather of tubes than of lines. 
However, this is not important because we are interested 
now in comparatively high frequencies u> of the order of 
a typical period of the orbit. Very long orbits would con- 
tribute to the function I(u>), Eq. (|6.3| ), at much smaller 
frequencies. The question about the finite thickness of 
the lines can aris e for times larger than the Ehrenfest 
time tE, Eq. (6.1), which we do not consider here. 

So, let us assume that we have got a periodic orbit. All 
supermatrices T and Q in Eqs. (16.91) are assumed to be 



on this orbit and we can write these equations in a more 
convenient form choosing the thickness of the paths and 
changing to the continuous limit along them. Taking the 
continuous limit along the paths the elementary length 
/ will disappear. The only quantity to be chosen is the 
elementary phase volume [Af2] in Eq. (6.6). 

In the d-dimensional space the density of states v can 
be written as 



dn 

de F 



Cl d d 



(2irh) a v F 



(6.10) 



where fid is the volume of the unit d-dimensional sphere. 

Then, using the relation between the surface of the 
unit sphere and its volume Sd = fldd we reduce the coef- 
ficient u\Q ] entering the free energy functional -F±[Q] in 
Eqs. to the form 



d-l 



(6.11) 



Now we have to choose the product n»=i (Apj_ Ar^_) an d 
we do this using a standard quasiclassical rule according 
to which we write 



(6.12) 



With Eq. fl6.12p we obtain 

u[AQ] = (2Trhv F y 1 



(6.13) 



Putt ing as everywhere before li = lwe can rewrite Eqs. 
.9|) as a sum over periodic orbits 



£/<(Q 44 (*)-l) 



{2w F y ^ 

x (Q 88 (x') + l))idxdx' 



(6.14) 



where < ... >i stands for the functional integral 



(->l 



J (...)exp(-Fi[Q])£)Q (6.15) 



and A = (lyV)^ 1 is the mean level spacing for the billiard 
under consideration. The one-dimensional free energy 
functional F\ [Q] for an orbit takes the form 



Fi [Q] = \str J 



AT (x) 



dT (x) i (lu + iS) 



dx 



2v F 



AQj dx 
(6.16) 



and is th e result of the adding of th e free energies F± [Q±] 
Eq.( |6.9| ). The sum over p in Eq. (6.14) means the sum 
over all periodic orbits and, in principle, Q (x) should de- 
pend on p. In order to simplify notations we omit writing 
this dependence explicitly. The integrals over x and x' 
are taken along the orbits (a certain direction is implied 
to be already chosen). 

The overall coefficient in the function al F\ [Q] is de- 
termined by the quasiclassical r ule, E q. ( |6.12| ). It is not 
difficult to understand that Eq. ( |6.12] ) is the only reason- 
able choice for the "thickness" of the classical paths. The 
functional F\ [Q] is multivalued because any replacement 
of the type T (x) — * T (x) h (x), where h (x) h (x) = 1 and 
[h (x) , A] = 0, does not change the supermatrix Q (x) but 
changes the functional F% [Q] as 



F 1 [Q] -> Fr [Q] 



(6.17) 



1 f „ f T , . dhi (x) T , . dh 2 (x) . . 
+ - / Str U a (x) ^ ; - h 2 (x) ^ ' ) d.r 



where h± and h 2 are the upper and lower diagonal blocks, 
respectively. Writing h m , m = 1, 2 as 



hm (x) = ( ex p(^K) ° 

V exp(ix m (x)) 



(6.18) 



<t>1n (*) 

-fa (x) J ' 



Xrn ( X ) 
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we write the change Aiq [Q] of the functional F\ [Q] as 



(6.19) 

where A<^ and are changes of the phases when 

going around the periodic orbit. These changes must 
be integer multiple of 2tt. The phases (x) and 4>m i x ) 
(as well as xtn ( x ) an d X m ( x )) are n °t independent of 
each other. If 0+ (x) = (2irk m /L p )x, then <f>~ n (x) = 
— (27rfc m /Lp) a; (fc m is an integer and L p is the length of 
a p-orbit). In a general case we can write 

<t>t(x)=±^x<f>t(x), ^(0) = ^(i P ) (6.20) 

and the same for Then, we obtain for the changes 
of the phases A0+ = -A<f>~, A x + = -A*", which 
reduces Eq. ( |6.1S ) to the form 



AFi [Q] = 2niM 



(6.21) 



where M is integer. 

With Eq. ( |6.21 ) we come to the result that, although 
the functional F\ [Q] is multivalued, the partition func- 
tion exp(— F\ [Q]) is not and one can integrate over Q 
in a standard way. In principle, we could proceed in the 
opposite way and determine the thickness of the classi- 
cal paths by demanding the partition function be single 
valued. Then, we would obtain the quasiclassical quan- 
tization rule, Eq. ( 6.12| ), automatically. 

Now , let us calculate the functional integral, Eq. 
( |6.14 ), using two different approaches. First, we compute 
the functional integral following the scheme of Ref .E3. We 
can us e diff eren t par ametrizations like those speci fied b y 
Eqs. ( |3.14| ) or (|4.28D . The parametrization, Eq. fl4.28| ), 
is very convenient because the Jacobian is equal to unity. 
At the same time, the ballistic cr-model contains also non- 
harmonic in P terms. In the quadratic approximation, 
one writes the free energy as 



F^[Q] = -Str I (AP(x) 



dP (x) i (uj + id) 



dx 



Vf 



P 1 (x)j dx 
(6.22) 



The supermatrices Q (x) entering Eq. ( |6. 14 ) should also 
be expanded up to quadratic terms. Calculating gaussian 
integrals we obtain for I (w) 



I (to) 



where 



2, ,2 



TT*V 



E E 

p m= — oo 



,0) 



2nm 



uj + iS 
vf 



, (6.23) 



are momenta corresponding to a p -periodic orbit, L p is 
the len gth of the orbit. 

Eq. (6.23) corresponds to the "perturbative part" of 
the level-level correlation function of Ref.E3, although 
here we sum over momenta on periodic orbits instead of 
summation over eigenvalues of the Perron- Frobenius op- 
erator. However, it is clear that Eq. (6.23) cannot corre- 
spond to a good perturbation theory because this expres- 
sion contains resonances at arbitrarily high frequencies. 

One might guess that next orders of the expansion in P 
had to be taken into account and the expansion would not 
be good at the resonances. Curiously enough, it is not 
so. One c an see immediately that the quadratic form , 
Eq. ( |6~22|) is exact in the parametrization, Eq. ( |3.14| ) 
(the overall coefficient is 2 times larger). At the same 
time, the part Q" o f the supermatrix Q commuting with 
A and entering Eq. (|6.14|) is exactly QH — 1 + 2P 2 . There- 
f ore, w e do not obtain any perturbative corrections to Eq. 
( |6.23|) (at the same time, the co ntribu tion of the Jaco- 
bian in the parametri zatio n, Eq. (3.14), is not as clear). 

Never thele ss, Eq. (|3.23| ) is not exact. The matrices P 
in Eq. (3.14) vary on non-trivial manifolds and extend- 
ing the integration over these matrices from — oo to +oo 
as it is implied in any gaussian integration is unjustif ied . 
To m ake the discussion simpler, let us rewrite Eqs. 



6.23) using the Poisson summation formula as 



T„A 



r(u) = i + e ( - 1 — Re E nexp ( w + iS ) T P n ) 



71=1 



(6.24) 



where T p = L^/v p is the period of the motion on the p - 
orbit. Eq. (6.24) corresponds to an expansion in pe riodic 
orbits of a classical flowEJ. Strictly speaking, Eq. ( 6.24 ) 
is different from what one writes for classical flows by 
absence of a factor containing the monodromy matrix. 
It is clear that in our simple consideration the stabil- 
ity of the periodic orbit is not taken into consideration. 
As any periodic orbit we consider has a finite thickness, 
the monodromy matrix would appear in a more accurate 
calculation. However, the aim of this chapter is only to 
clarify the origin of the repetition problem and therefore 
we use the simplest approximation. 

The factor n in front of the exponential is a charac- 
teristic feature off-expansions for classical flows (see e.g. 
Eq. (52) of Ref.c2l which leads to this dependence af- 
ter taking the logarithm of both parts and taking second 
derivative in s). In other words, we have now an expan- 
sion in periodic orbits of the Perron-Frobpaius operator 
and this corresponds to the result of Ref.EHl in the limit 
of the vanishing regularizer. 

However, although the perturbative approximation of 
Ref.E3 works very well in the diffusion limit, we do not 
see any justification for it in the ballistic limit. There- 
fore, we should try to calculate the integral in Eqs. (3.14, 
6.15) without using this approximation 



Fortunately, the functional integral in Eq. ( 6.14 , 6.15| ) 
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can be calculated exactly even easier than approximately. 
The free energy functional Pi [Q] entering these equations 
corresponds to a one-dimensional ring without any impu- 
rities, provided the averaging over the spectrum has been 
performed. This energy averaging is neces sary to get a 
smooth quasiclassical function g n (r), Eq. ( 2.20 ). With- 
out the averaging this function would not be smooth due 
to quantization of the energy levels in the ring. So, we 
conclude that the calculation of functional integral over 
Q with the free energy functional Pi [Q] is equivalent to 
calculation of the averaged level-level correlation function 
R p (uj) for a clean electron system on a ring. 

The level-level correlation function Ri p (oj) for such a 
ring can be written as 

oo 

Rip (uj) = (viL p y 2 ^2 (S (e - w - e (q m )) 

m,m' — — oo 

x5(£-e{q m ,))) e (6.25) 

where (. . .) E means the averaging over e, q rn = 2nm/ L p 
and vi — (ttvf) 1 is one-dimensional density of states. 
The spectrum e (q) can be, as usual, linearized 



e{q m ) = q -kr^~v F (\q m \-p* 
2m \ 



(p) 



(6.26) 



Using the Poisson formula and Eq. ( 6.26 ) we write 

D j- OO 

6 (e - e (q m )) ^ ex P (* ( 27 "4 P) + bTp) 

(6.27) 



m— 



n— — oo 



where n 



(/-<) 



Substituting Eq. (6.27) into Eq. (6.25) we obtain 



oo 

R lp (u))= Y cxp (i(2mi ( p + eT p )(n - n') + ujT p n'^ 



n.n — — OO 



(6.28) 

After averaging over e, o nly t he term with n = n' gives 
the contribution in Eq. (6.28). On the other hand, re- 
peating all the steps of th e derivation of the cr-model, we 
come to Eqs. (3.14, 6.15) with the only difference that 
we should replace the mean level spacing A of the en- 
tire system by the level spacing of the one-dimensional 
ring A p = ir/T p . This allows us to write the level-level 
correlation function R (uj) of the quantum billiard under 
consideration as 



i?M = i + ]T 



Re exp (i (uj + id) T p n) 



(6.29) 



Comparing Eqs. (3.24) and ( 3.29| ) with each other we see 
that the only difference between them is the presence of 



the prefactor n in Eq. (|J3). So, we conclude that the 
factor n in the expansion in periodic orbits i*. a conse- 
quence of a unjustified approximation of Ref.O and this 
solves the problem of repetitionsE3. The assumption of 
Rcf.tL3 that the regularizer can be put to zero at the end 
of calculationSj-does not seem to be correct. We believe 
(following RefO) that the presence of a finite regularizer 
is inevitable in a quantum system and a very important 
problem is to calculate it. 

The problem of repetitions was discussed recently for 
weak scatterers in RefO where the problem was related 
to the question of a possibility of separating 4 point cor- 
relation functions into two diffusons. From the above 
discussion, we see that the problem is even more deli- 
cate because Eq. ( 6.23 ) is perturbatively exact and the 
difference comes from oscillating exponentials. 

In order to understand better what has bee n ne- 
glected in our derivation, we compare Eq. ( |6.29 ) with 
a corresponding diagonal contribution obtained from the 
Gutzwiller trace formula, (see, e.g.Eil) 



RW ( a 



T 2 

1 H 



EE*- 

p n—l 



T' 2 
p 



p 



{ 2TmT p 



(6.30) 



where g p is the action-multiplicity of the p primitive 
orbit (for the orthogonal ensemble g p — 2), and M p is 
the monodromy matrix that describes the flow linearized 
in its vicinity. In Eq. (6.30), Tjj = 2ir/A is the Heisen- 
berg time, T p is the orbit period, and x = ui/A. Except 

for the factor \M£ - if 1 Eqs. ( |6~29| ) and |T3C| ) agree. 
As we carried out computation without any regularizer 
like the one of RefO, we conclude again that its pres- 
ence is absolutely necessary and it must be related to the 
monodromy matrix M p . We see that without this term 
correlations between orbits do not exist and one cannot 
pass to the universal limit when lowering the frequency uj. 
In the language of the field theoretical approach, one can- 
not reduce the ballistic cr-model to the zero-dimensional 
one without this regularizing term that must describe 
quantum diffraction on irregularities of the boundary. 
Adding a term like 



F reg [Q] - Str J n (r) drdn (6.31 



where /3 n (r) is a function in the phase space, we may 
obtain effectively a coupling between the periodic orbits 
of the type 

-^Ppp'Str J ( Xp ) Q( p ') (ay) dx p dx p , (6.32) 
p,p> 

which resembles coupling between grains in a granular 
system. For small j3 pp i the orbits are not coupled and 
we have separate periodic orbits. As the coupling (3 PtP i 
grows, the relative fluctuations of Q&> with respect to 
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each other get suppressed and one needs to consider ro- 
tations of the system as the whole. Then, one obtains the 
zero-dimensional er-model and, hence, the Wigner-Dyson 
statistics. This is a scenario for a granular metal and 
we believe that it is relevant for the quantum billiard, 
the role of the grains being played by separate periodic 
orbits. 



VII. DISCUSSION. 

In the present work, we made an attempt to put the 
field theoretical approach to systems with a long range 
disorder on a solid basis. The conventional method of 
derivation of the supermatrix er-modelEI is based on sin- 
gling out slow modes, performing Hubbard-Stratonovich 
transformation and using a saddle-point approximation. 
Although this approach worked well for a short range 
disorder, its validity is not justified for a long range dis- 
order and quantum billiards. We suggested a scheme 
that allows us to overcome these difficulties an d de rive a 
modified non- linear ballistic er-model (see Eq.(L13)). 



The method resembles the approach of RefJ13 and is 
based on writing quasiclassical equations for generalized 
Green functions. At the same time, the quasiclassical 
equations are written for non-averaged over the long 
range potential quantities and singling out slow modes is 
performed only for a part originating from a short range 
disorder. In addition, the short range disorder does not 
play an important role and can be put to zero. The cru- 
cial step of the derivation is that the solution of the qua- 
siclassical equations can be found exactly, which is a con- 
sequence of the supersymmetric structure of the Green 
functions. This possibility was overlooked in the previous 
studyH. 

The scheme developed now leads to a considerable 
progress in describing disordered systems with long range 
disorder because the derivation is applicable for all 
lengths exceeding the wavelength Xp. Representing the 
solution of the quasiclassical equations and also the parti- 
tion function for the electron Lagrangian with sources in 
terms of a functional integral over supermatrices Q n (r), 
n 2 = 1, with the constraint Q\ (r) = 1 we were able to 
average over the disorder exactly and obtain a ballistic 
er-model in a new form that has.not been written before. 
The so called "mode locking" Ej problem does not arise 
here because the eigenvalues of the supermatrix Q n (r) 
are fixed by the construction and do not fluctuate. The 
method suggested resembles the method of bosonization, 
well known in field theory, see e.g. a booko, when a 
fcrmionic system is replaced by a bosonic one. In our ap- 
proach, we also replace the electron system by a system 
of ballistic excitations that can be considered as quasi- 
particles. At large scales, these quasiparticles are well 
known diffusons and cooperons. In analogy, our scheme 
can be called superbosonization. 

For weak scatterers, there should exist 2 more scales: 



the Lapunov length I 1 and the transport mean free path 
Itr- The single particle mean free path I does not ap- 
pear in our consideration (actually, we do not consider 
one-particle Green functions at different points restrict- 
ing our study to gauge invariant quantities). Integrating 
out degrees of freedom related to distances smaller than 
the Lapunov length I 1 we obtained a reduced ballistic er- 
model. A propagator describing small fluctuations within 
this reduced er-model corresponds to the kinetic Boltz- 
mann equation with a collision term. Integrating further 
on scales up to the transport mean free path l tr we obtain 
the standard diffusive er-model. 

Trying different calculational schemes we conclude that 
one can do perturbative calculations with the ballistic er- 
model only at scales exceeding the Lapunov length (we 
call this range—j'collision region"). At smaller lengths 
(following RefJiZI we call this range Lapunov region) no 
perturbation expansions in diffusons and cooperons are 
possible. In this region one can carry out calculations 
deriving equations for correlation functions and investi- 
gating them in different approximations. 

It seems that an infinite system with a weak long range 
disorder is adequately described by the ballistic er-model 
we have derived. However, when describing quantum bil- 
liards, our approach is not accurate near the boundaries, 
where the quasiclassical approximation may not be used 
(turning points). We believe that a more accurate deriva- 
tion may result in a new term in the er-mod&L|-This term 
was suggested phenomenologically in Refs.EZlEl but has 
not been derived yet microscopically. Its presence seems 
to be absolutely necessary because it must introduce a 
new scale: the Ehrenfest time tp- One may not put 
the regularizer to zero at the end of calculations. We 
have demonstrated that neglecting such a term we re- 
duced the ballistic er-model for the billiard to ballistic 
er- models for periodic orbits. Proceeding in this way we 
demonstrated explicitly where the contradiction between 
the workli3 and Refill (repetition problem) comes from. 
Our conclusion is that the representation of the level- 
level correlation function in terms of eigenvalues of the 
Perron- Frobenius operator suggested in Ref.Ej is not jus- 
tified in the ballistic case. This approach is valid only if 
there are no repetitions but this would rather correspond 
the diffusive case. In the opposite limit, one comes to 
a description in terms of periodic orbits without corre- 
lations between actions of different psbits. This is the 
region where the description of Ref.ta may be applica- 
ble. At the same time, the Ehrenfest time can hardly be 
identified on the basis of the trace formula and. therefore 
the limits of applicability of the result of Rcf£3 have not 
been specified. The hypothetical regularizer seems to be 
related to the monodromy matrix entering the Gutzwiller 
trace formula. 

We believe that the field theoretical approach pre- 
sented here and the formalism based on the Gutzwiller 
trace formula can be complementary to each other de- 
scribing quantum systems in different regions of parame- 
ters. At times smaller than the Ehrenfest time, the trace 
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formula can be more convenient. However, at a larger 
time, trying to extract physical quantities from the trace 
formula does not make much sense because the cr-model is 
a much more convenient tool for such calculations. This 
concerns especially the universal limit where the cr-model 
approach leads for most correlation functions to definite 
integrals that can be computed rather easily. 

It is important to notice, that .diagrammatic expan- 
sions like those attempted p. Ref.EJ can hardly be suc- 
cessful. The authors of RefO found that the results de- 
pended crucially on the way how the ultraviolet cutoff 
was introduced. Now we understand that the ultraviolet 
cutoff must be imposed by the requirement that the in- 
tegration is performed over the manifold Q\ (r) = 1 in 
an invariant way. Any artificial ultraviolet cutoffs in the 
perturbation theory would correspond to a violation the 
rotational invariance in the space of the supermatrices 
Q n (r) and lead to wrong results. At the same time, it 
is not clear how to develop a perturbation theory in an 
invariant way. 

An important question aL aa i-aMer aging procedure was 
discussed in several workslll'EJO'ta. In all these publi- 
cations an opinion was expressed that an averaging over 
energy was not sufficient for study of quantum chaos in 
quantum systems and different types of an additional 
averaging were suggested. We do not agree with this 
point of view because, in our derivation of the quasiclas- 
sical equations, averaging over the energy allowed us to 
smooth generalized Green functions and this was all we 
needed. The only condition is that the averaging should 
be performed in an interval of energies much exceeding 
the mean level spacing. 

The souro&,©£||the discrepancy is simple: the authors 
of the worksEj'tSlla used the saddle point approximation 
and the expansion in gradients. Although rtfae saddle 
point approximation was not necessary in Ref.t2l, the gra- 
dient expansion still had to be carried out. Therefore, an 
additional averaging was necessary to justify these ap- 
proximation. Since we do not do such approximations, 
no additional averaging is needed in our scheme. 

In conclusion, our approach enables us to carry out 
calculations for long range disorder and ch aos in a reli- 
able way. Stil l, a derivation of a new term (6.31) in the 



cr-model (3.13) describing quantum diffraction on bound- 
aries has to be done to make the theory complete but we 
believe that this is not impossible. 
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APPENDIX A: BOUNDARY CONDITIONS 

In this Appendix, we derive boundary conditions for 
the boundary of the sample. We describe the boundary 
by an external field u_e(r) which is negligible inside the 
sample and grows sharply at the surface of the sample. 
In such a situation, the electron wave function decays 
fast outside the sample and, in the limit of infinite po- 
tential walls, one can just put the wave function equal to 
zero at the boundary. Unfortunately, such a boundary 
condition is not very helpful because in the quasiclassical 
approximation we use it looses its validity at a distance 
of several wave lengths from the boundary and matching 
the wave functions in the bulk and at the boundary is 
necessary. 

In order to find effective boundary conditions for qua- 
siclassical Green functions we follow methods well devel- 
oped in superconductivity theoryEI. First, we write the 
supermatrix G"^(r;r') in a form of a sum over eigensu- 
pervectors Vfc( r ) 



G af3 (i 



k 



Eh- £ 



(Al) 



satisfying the Schrodinger equation 



, s . u + i5 . . 
H 0r + u (r) + A — ViJ{r) 



i>k (r) = e k ip k (r) 

(A2) 



We assume that the potential u (r) in Eq. ( |A2j ) contains 
not only the impurity field but also the potential ub (r) 
describing the boundary. The conjugated equation can 
be written as 



^fc(r) 



/ N .Ul+iS ' T f \ 

H Qr + u(r) + A — hiJ(r) 



s k Mr) (A3) 



The summation in Eq. QAip should be performed 
over the complete set of eigenfuctions, so that 
12k tPk MV'f ( r ') = 8 af3 5(r — r'). It means that the choice 
of the set of eigenf unctions ip^ (r) and the operation of 
the conjugation must conform with each other. For ex- 
ample, at distances from the boundary much larger than 
the wavelength Xp but much smaller than the radius b 
of the random potential, we choose the cigenfunctions 
ip k (r) in a form of plane waves 



^p(r) 



(A4) 



where ip p is a normalized vector from the superspace. 
Then, we have to adjust the definition of the conjugation 
written in the bookO and add to it the momentum inver- 
sion p — » —p. If the motion cannot be treated as free 
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changing the sign of the momentum has to be generalized 
by replacing this op era tion by the time reversal. Since 
the energy Ek in Eq. ( [A2] ) remains the same after the time 
reversal, the spectra l exp ansion written in Eq.(Al) is in 
agreement with Eq.(2.17). Below, we will use this rela- 
tion between supervectors tp a (r) and those conjugated to 
them. 

Now we introduce local coordinates (z, r M ) in the vicin- 
ity of the boundary. The coordinate r M is a coordinate 
along the boundary surface and z- is a distance between 
a given point and the surface. Points on the surface have 
the coordinates (0,r M ). If the boundary is rough, the co- 
ordinate system (z,r u ) is not very useful. However, if the 
boundary is smooth, which means that the derivative of 
the field fig(r) along it is small in comparison to that 
in the perpendicular direction, the coordinates (z, r„) are 
very convenient for the quasiclassical approximation. 

If the radius of the curvature of the boundary is large 
the electron wave function in the vicinity of the boundary 
can be represented as a sum of one-dimensional solutions 
with respect to the z-direction with amplitudes slowly 
dependent on the coordinates r M . Since we need to know 
the wave function in the domain in which the potential 
ub{y) vanishes, we can write the asymptotic form of the 
wave functions as 

^ fe (r) = e*P" r " (e*-Vfc,c(r„) + e^Vfe.r (r„)) (A5) 

where tpk,c{ r ") and <^fc, r (rn) are some slowly varying su- 
pervector functions depending on the potential ub (r) 
that play the role of amplitudes of the coming and re- 
flecting waves respectively. They comprise the minimal 
knowledge about the potential us(r) that is needed to 
find the required boundary conditions. 

If the boundary is impenetrable and the potential re- 
flects all waves, then the amplitudes tfk,c{ r ")' <Pfc,r( r n) 
can be determined from the condition that the compo- 
nent of the current perpendicular to it is equal to zero. 
This condition is valid not only in the region in which the 
potential ug(r) is relevant but also in the quasiclassical 
region because the current cannot considerably change at 
distances of the order of the wavelength. 

An expression for the current can be obtained in a 
standard way fro m t he par ticle conservation law that 
follows from Eqs.(A2), (A3). Its z-component perpen- 
dicular to the surface is proport iona l to the difference 
dz'tpip ~ i^9 z ip. Substituting Eq. JA5|) and its conjugate 
into the difference and putting the result to zero we find 

ipk,c( r ») l Pk,c( r ») = 9k,r( r ») l Pk,r( r »)- A relation between 
the amplitudes of the coming and reflected waves in the 
case of the impenetrable boundary can also be estab- 
lished by the demand that they should transform one 
into the other by the time reversal. A general expression 
that satisfies it can be chosen in the form: 

V fc (r) = e*P" r " ( e <P.<*— - e -*P.(»-*o) ^( r „) (A6) 

where zq determines an unknown phase tha t ca n be found 
only by matching the function t[>k (r), Eq. (A6), with the 



corresponding decaying asymptotics at the opposite side 
of the turning point (its value is of the order Af). This is 
generally not an easy task but, fortunately, the parameter 
zq is not important for finding the boundary conditions 
for the Green functions. 



S ubst ituting Eq.(A6) into the spectral expansion, 
Eq.(Al), we find an expression for the matrix G a/3 (r; r') 
at the boundary. To determine the matrix g n (v), Eq. 
fl2.20| ), it is necessary to carry out both the summation 
in £ (the system may be finite and we have to sum instead 
of integrating over £) and averaging over the energy. The 
latter is absolutely necessary because only this averaging 
guarantees vanishing of all terms containing the products 
ippn (r + r') in the exponents (n is a unit vector paral- 
lel to p). As soon as the terms containing ippn (r + r') 
vanish, the parameter zq drops out. Then, the Green 
function G Q ^(r;r') can be written in the vicinity of the 
boundary as 



(A7) 



xe ip„(r„-r:,)( e -ip,(*-*') + e +M z -*'))/ fe (r ll ,r' ll )) £ 

where f k (r„, rf,) = p k (r, and („.) E stands for av- 

eraging over the energy. 



Carrying out in Eq. (A7) summation over £ and aver- 
aging over the energy, which is equivalent to integration 
over £, and Fourier transforming with respect to r — r' 
we obtain the quasiclassi cal function <?„ (r) . The func- 
tion obtained from Eq. (A7) is a slow function of the 
coordinate r + r' and does not change at the boundary 
under the replacement — > — n±. Hence, we come to 
the boundary condition for the quasiclassical Green func- 
tions g n (r) 



r = (z = 0,r„) 



(A8) 



where is the component of the vec tor n perpendicu- 
lar to the surface. Although Eq. ( |A8|) is rather simple, 
it has not been written in previous works on the ballistic 
tr-models. 
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